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ABSTRACT 

We discuss structural and kinematical properties of the stellar halo and the old 
globular cluster system (GCS) in the Large Magellanic Cloud (LMC) based on numer- 
ical simulations of the LMC formation. We particularly discuss the observed possible 
GCS's rotational kinematics {V/a ~ 2) that appears to be significantly different from 
the stellar halo's one with a large velocity dispersion 50 km s~^). We consider that 
both halo field stars and old GCs can originate from low-mass subhalos virialized at 
high redshifts {z > 6). We investigate the final dynamical properties of the two old 
components in the LMC's halo formed from merging of low-mass subhalos with field 
stars and GCs. We find that the GCS composed of old globular clusters (GCs) formed 
at high redshifts {z > 6) has little rotation {V/a ~ 0.4) and structure and kinematics 
similar to those of the stellar halo. This inconsistency between the simulated GCS's 
kinematics and the observed one is found to be seen in models with different parame- 
ters. This inconsistency therefore implies that if old, metal-poor GCs in the LMC have 
rotational kinematics, they are highly unlikely to originate from the low-mass subha- 
los that formed the stellar halo. We thus discuss a scenario in which the stellar halo 
was formed from low-mass subhalos with no/few GCs whereas the GCS was formed 
at the very early epoch of the LMC's disk formation via dissipative minor and major 
merging of gas-rich subhalos and gas infall. We also discuss whether old GCs in the 
LMC can be slightly younger than the Galactic counterparts. We suggest that there 
can be a threshold subhalo mass above which GCs can be formed within subhalos at 
high redshifts and thus that this threshold causes differences in physical properties 
between stellar halos and GCSs in less luminous galaxies like the LMC. 
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1 INTRODUCTION 

Structural, kinematical, chemical properties of the Large 
Magellanic Cloud (LMC) have been investigated by many 
authors concerning different stellar populations and gaseous 
components and suggested to have valuable information on 
formation and evolution history of the LMC (Hartwick & 
Cowley 1988; Meatheringham et al. 1988; Irwin 1991; Luks 
& Rohlfs 1992; Kunkel et al. 1997; Graff et al. 2000; Olsen 
& Salyk 2002; van den Marel et al. 2002; Cioni & Habing 
2003; Staveley-smith et al. 2003; Cole et al. 2005). Previous 
numerical studies on dynamical and chemical evolution of 
the LMC compared the results with these observations and 
thereby tried to provide reasonable physical explanations 
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for the observations (e.g., Bekki et al. 2004; Bekki & Chiba 
2005). Since these studies focused mainly on the origin of 
disk components of the LMC, the origin of the old stellar 
halo and the globular cluster system (GCS) composed of 
old, metal-poor globular clusters (GCs) remains unclear. 

Recent observations have reported that the stellar halo 
of the LMC has a projected radial density profile similar 
to an exponential one (e.g., Alves 2004) and a larger ve- 
locity dispersion with possibly little rotation (Minniti et al. 
2003; Gratton et al 2004). The dynamically hot nature of 
the LMC's stellar halo has been confirmed by Borissova et 
al. (2006), which found that the stellar halo composed of RR 
Lyrae stars has the mean velocity dispersion of ~ 50 ± 2 km 
s~^ and a Gaussian metallicity distributions function with 
mean [Fe/H] = -1.53 ± 0.02 dex. Subramaniam (2006) in- 
vestigated spatial distributions of RR Lyrae stars from the 
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Table 1. Model parameters and a brief summary of the results 



model 


Mt (x1O1OM0)_^ 






N ■ 

mm ^ 


^trun ^ 


V <T ''fS — 


(Z) 9 
y o- -'gc — 


Standard 


6.0 


0.08 


0.39 


32 


15 


0.34 


0.39 


Low density 


6.0 


0.08 


0.19 


32 


15 


0.60 


0.56 


High threshold 


6.0 


0.08 


0.39 


1000 


15 


0.28 


0.22 


Low-z truncation 


6.0 


0.08 


0.39 


32 


10 


0.19 


0.21 



Initial total masses. 
^ Initial spin parameters. 

Initial over densities. 

Minimum mass for halo identification. 
^ The epoch of truncation of GC formation. 

f Final — for field stars (FS), where V and cr are the maximum rotational velocity and the central velocity dispersion, respectively in 
simulated LMC halo. The man value averaged for three projections {x-y, x-z, and y-z) is shown. 
9 y- for GCs. 



catalogue by Soszynski et al. (2003) and found that the inner 
stellar halo can have a disky density distribution. 

Previous observational studies on the kinematics of the 
GCS composed of old, metal-poor GCs in the LMC sug- 
gested that the GCS has a disky distribution and rotational 
kinematics with a small velocity dispersion (e.g.. Freeman et 
al. 1983; Kinman et al. 1991; Schommer et al. 1992; Grochol- 
ski et al. 2006). It is suggested that it could be difficult to 
make a robust conclusion on the presence of the rotational 
kinematics owing to the small number (13) of the GCs (van 
den Bergh 2004). The total number of old, metal-poor GCs 
with [Fe/H] < —1.3 and well determined radial velocities 
are five in the sample of GCs by Grocholski et al. (2006), 
which implies that old GCs can have different kinematics 
from young and intermediate-age GCs with rotational kine- 
matics. The observed possible rotational kinematics with 
V/a ~ 2 (Schommer et al. 1992) however suggests that there 
is a significant kinematical difference in the stellar halo and 
the GCS composed of old GCs in the LMC. No theoretical 
explanations however have been proposed for the origin of 
this kinematical difference. 

The purpose of this paper is to investigate the origin of 
the observed kinematical difference between the stellar halo 
and the GCS composed of old, metal-poor GCs in the LMC 
based on numerical simulations of the LMC formation. We 
here consider that the stellar halo and the GCS are formed 
from hierarchical merging of subhalos that are virialized be- 
fore reionization and have both field stars and GCs. We 
therefore perform numerical simulations for the formation 
of the LMC based on the cold dark matter (CDM) model 
and thereby investigate structure and kinematics of the sim- 
ulated stellar halo and GCS. By comparing the simulated 
kinematical properties of the stellar halo and the GCS with 
observations, we try to discuss the relationship between the 
stellar halo formation and GC one in the early history of the 
LMC. 

The plan of the paper is as follows: In the next section, 
we describe our numerical models for the stellar halo and 
GC formation of the LMC. In §3, we present the numeri- 
cal results mainly on the kinematics of the simulated halo 
and GCS for variously different models. In §4, we discuss a 
promising scenario explaining the observe kinematical differ- 
ence between the old stellar halo and the GCS in the LMC 
We summarize our conclusions in §5. 



2 THE MODEL 
2.1 GC formation 

This paper is the first step toward better understand- 
ing structural, kinematical, and chemical properties of the 
LMC's GCS in a comprehensive way. Therefore we adopt 
a more idealized model of GC formation in the LMC that 
is assumed to be formed from hierarchical merging of low- 
mass subhalos. We adopt a model in which GCs (and field 
stars) are formed within subhalos virialized at high redshifts 
{z > 6). Although the adopted model of GC formation 
within subhalos at high redshifts has not been confirmed 
observationally (Brodie & Strader 2006 for different models 
of GC formation), recent numerical simulations based on the 
model have successfully explained some fundamental obser- 
vations such as radial density profiles of GCSs (e.g., Santos 
2003; Bekki 2005). 

We consider that star formation and thus GC one can 
proceed at high redshifts only in virialized dark matter halos 
before reionization. The physical reason for the suppression 
of star formation by reionization is that ultraviolent back- 
ground radiation in a reionized universe can significantly 
reduce the total amount of cold HI gas and molecular one 
(through photoevaporation/photoionization of the gas) that 
are observed to be indispensable for active star formation in 
low-mass galaxies (e.g., Young and Lo 1997). Recent high- 
resolution simulations on this issue (Susa & Umemura, 2004) 
have confirmed that significant suppression of the formation 
of cold gas can lead to the suppression of star formation in 
dwarf galaxies embedded in dark matter halos, in particu- 
lar, lower-mass dwarfs. We here focus exclusively on very 
old, metal-poor GCs in the LMC, which could have formed 
at very high redshifts: it is, however, observationally unclear 
whether the GCs were formed before reionization. Therefore, 
it is reasonable to adopt the above assumption that GC for- 
mation can proceed in dark matter subhalos virialized before 
reionization. 

Massive and dense star clusters like GCs are suggested 
to form in GMCs surrounded by high-pressure interstellar 
medium (ISM) in galaxies (Elmegreen & Efremov 1997) . We 
therefore assume that GCs are formed in the very central re- 
gions of subhalos, where gaseous pressure should be so high 
owing to their deep gravitational potentials. During destruc- 
tion of subhalos with GCs in the hierarchical growth of the 
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LMC, the GCs are stripped and dispersed into the halo to 
become the halo GCs in the LMC. This formation processes 
of GCs from low-mass galactic systems was demonstrated by 
previous numerical simulations (e.g., Bekki & EYeeman 2003; 
Mizutani et al. 2003; Bekki & Chiba 2004), which suggests 
that the adopted assumption is quite reasonable. 

We perform purely dissipationless simulations on 
galaxy-scale halo formation via hierarchical merging of sub- 
halos with field stars and GCs. In these dissipationless simu- 
lations, we first identify possible formation sites of field stars 
and GCs in low-mass halos at high z and then follow their 
evolution during hierarchical merging of the halos till z = 0. 
The final structural and kincmatical properties of the stellar 
halo and GCS in the simulated LMC arc thus determined by 
the details of merging histories of subhalos with field stars 
and GCs. We consider that as long as the formation sites of 
field stars and GCs in low-mass halos are properly modeled, 
the present dissipationless models allow us to derive physical 
properties of the stellar halo and the GCS in the LMC in a 
reasonable way. The present simulations arc different from 
our previous chemodynamical ones (Bekki & Chiba 2000; 
2001) which investigated both dynamical and chemical evo- 
lution of forming galaxies in order to reproduce structures, 
kinematics, and chemical properties for the Galactic field 
stars in a self-consistent manner. We consider that it is cur- 
rently difficult and numerically costly to construct the fully 
self-consistent chemodynamical models in which formation 
sites both for field stars and GCs can be directly derived for 
gaseous regions of low-mass halos. Accordingly, the present 
study is the first step for better modeling star and GC for- 
mation in the LMC. Fully self-consistent chemodynamical 
simulations with reasonable and realistic models for star and 
GC formation from GMCs will be done in our future works. 

2.2 Identification of field star and GC particles in 
hierarchical galaxy formation 

We simulate the formation of galaxy-scale halos in a ACDM 
Universe with Q, = 0.3, A = 0.7, Ho = 70 km s"^ Mpc"\ 
and (78 = 0.9, and thereby investigate merging/accretion 
histories of subhalos that can contain low mass dwarfs. The 
way to set up initial conditions for the numerical simula- 
tions is essentially the same as that adopted by Katz & 
Gunn (1991) and Steinmetz & MuUer (1995). We consider 
an isolated homogeneous, rigidly rotating sphere, on which 
small-scale fluctuations according to a CDM power spec- 
trum are superimposed. The initial total mass (iWt), radius, 
initial overdensity {5i), and spin parameter (A) are set to be 
free parameters. 

Although we investigate variously different models with 
6 X IO'^Mq < Mt «; 6.0 X 1O"M0, A = 0.08, 0.12, and 0.16, 
and Si = 0.19 and 0.39, we mainly show the results of the 
"LMC" models with Mt = 6.0 x 10^°Mo and A = 0.08. 
The choice of A = 0.08 arc demonstrated to be quite reason- 
able for late-type disk galaxies in previous CDM simulations 
(e.g., Katz & Gurm 1991; Steinmetz & Miiller 1995; Bekki & 
Chiba 2000, 2001), and accordingly we consider that models 
with A = 0.08 are also reasonable for the present simula- 
tions for less luminous disk systems like the LMC. Si = 0.39 
is chosen such that the final central velocity dispersions of 
"stellar" components are similar to the observed ones. 

The low-mass (Mt = 6 x lO^M©) models show stellar 



halos with velocity dispersions being significantly lower than 
the observed one. The high-mass (Mt = 6 X 10" Mo) models 
are investigated in order that we can compare the results of 
the LMC model with those of "the Galaxy" one and thereby 
discuss the differences in dynamical properties between the 
LMC's GCS and the Galactic GCS. These results wiU be dis- 
cussed in a wider context of stellar halo and GCS formation 
in galaxies with different Hubble types (Bekki & Chiba 2007, 
in preparation). The details of parameter values in the sim- 
ulations including those related to simulation methods (e.g., 
softening lengths) are described later. 

We start the coUisionless simulation at ^start (=30) and 
follow it till 2:cnd (=1) to identify virialized subhalos with 
the densities larger than 170pc(z), where Pc{z) is the crit- 
ical density of the universe, at a redshift z. This 170pc{z) 
corresponds to the mean mass density of a collapsed and 
thus gravitationally bound object at z (e.g., Padmanabhan 
1993). The minimum number of particles within a virialized 
subhalo (A'min) is set to be 32 corresponding to the mass res- 
olution of 3.8 X 10® Mq for the LMC models. This number 
of 32 is chosen so that we can find a virialized object at a 
given « in a robust manner. The mass resolution of 1.2 x 10^ 
AIq is chosen such that the masses of GC particles in the 
simulations can be consistent with the observed typical GC 
mass (~ 10^ Mq) in the Galaxy. 

For each individual virialized subhalo with the virialized 
redshift of Zvir, we estimate a radius (rb) within which 20 % 
of the total mass is included, and then the particles within 
rb are labeled as "baryonic" particles. This procedure for 
defining baryonic particles is based on the assumption that 
energy dissipation via radiative cooling allows baryon to fall 
into the deepest potential well of dark-matter halos. Such 
baryonic particles in a subhalo will be regarded as candi- 
date "stellar" particles to form stellax halos and GCSs in the 
later dynamical stage, if the subhalo is later destroyed and 
baryonic particles initially within the subhalo is dispersed 
into the galactic halo region. Thus, the present dissipation- 
less models track the formation of stellar halos and GCSs 
via hierarchical merging of subhalos, although the models 
are not adequate to the study of star formation histories in 
subhalos (as was done in our previous studies, e.g., Bekki & 
Chiba 2001). 

Stellar particles within rb in a subhalo are divided into 
field star ("FS") particles and "GC" ones accordingly to 
their locations with respect to the center of the subhalo. 
The stellar particle in the very center of a subhalo is identi- 
fied as GC particle whereas stellar particles other than the 
GC particle within rb are identified as FS ones. The initial 
distributions of FS particles are thus more diffuse than those 
of the GC particles in virialized subhalos. Massive, compact 
star clusters like GCs are suggested to be formed in extraor- 
dinary high-pressure regions, such as the centers of low-mass 
galaxies (e.g., Elmegreen 2004). Unbound or weakly bound 
star clusters, which can evolve into field stars after their dis- 
integration, can be formed in outer regions of galaxies where 
gas density and pressure arc low (e.g., Elmegreen 2004). We 
thus consider that the adopted model for the distributions 
of the two old stellar components is reasonable. The dif- 
ferences in initial positions with respect to the ccutcrs of 
subhalos between FS and GC particles can cause differences 
in the final distribution with respect to a galaxy-scale halo 
formed from the subhalos between these particles. 
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In the present model, each subhalo is assumed to have 
only one GC particle that is located in the center of the halo. 
One of the main reasons for this is that GCs in the present 
simulations are considered to be formed in nuclear regions 
of low-mass halos (i.e., GCs are initially either nuclear star 
clusters or stellar galactic nuclei) : this assumption is similar 
to the scenario proposed by Zinnecker et al. (1988). Further- 
more, our analytical arguments suggest that low-mass dark 
halos with the masses of 2 x IO^Mq, the luminous masses 
being 10% of the halos, and specific frequencies of GCs be- 
ing 5 can have one GC. Given that the halos identified as 
being virialized before reionization at each time step in the 
simulations mostly have masses less that 2 x 10® Mq, the 
above assumption of one GC in a halo can be reasonable. 

If the initial number of GC particles in each halo is 
increased, the initial difference in spatial distributions be- 
tween halo and GC particles becomes less remarkable in the 
halo: the final differences in structural and kinematical prop- 
erties between the stellar halo and the GCS in the LMC 
therefore becomes less remarkable. Thus the adopted mod- 
els with each subhalo having only one nuclear GC particle is 
regarded as those showing maximum possible differences in 
final dynamical properties between the stellar halo and the 
GCS in the present LMC model. 

2.3 Truncation of GC formation 

Previous theoretical studies have demonstrated that ultra- 
violet background radiation in a reionized universe can sig- 
nificantly reduce the total amount of cold HI and molecular 
gas that are observed to be indispensable for galactic ac- 
tive star formation (e.g., Susa & Umemura 2004). In order 
to investigate this suppression effects of star and GC for- 
mation on the final structural and kinematics properties of 
the simulated stellar halos and GCSs, we adopt the following 
idealized assumption: // a subhalo is virialized after the com- 
pletion of the reionization (zrdo-n), star and GC formation 
IS totally suppressed in such a subhalo. Then, hypothetical 
baryonic/stellar particles in the subhalos with z^h < Zroion 
will not be identified as FS or GC particles in the later stage, 
but those in the subhalos with Zvir ^ ^rcion will be regarded 
as progenitors of visible stellar halos. 

Recent WMAP ( Wilkinson Microwave Anisotropy 
Probe) observations have shown that plausible Zioion ranges 
from 11 to 30 (Spergel et al. 2003; Kogut et al. 2003) whereas 
quasar absorption line studies give the lower limit of 6.4 for 
Zreion (Fan et al. 2003). Guided by these observations, we in- 
vestigate the models with Zreion = (no reionization), 6, 10, 
15, and 20. The adopted picture of single epoch of reioniza- 
tion might well be somewhat oversimplified and less realistic, 
however, this idealized model can help us to elucidate some 
essential ingredients of the reionization effects on stellar halo 
and GC formation. For convenience, the epoch of the trun- 
cation of GC formation by reionization is denoted as ztrun 
in the following. 

We show the results of the models with 2:trun = 10 and 
15, firstly because these models explain structural proper- 
ties of the Galactic stellar halo and GCS (Bekki 2005; Bekki 
& Chiba 2005) and secondly because models with higher 
Ztrun can also better explain the observed properties of GCSs 
in early-type galaxies (Bekki et al. 2007). Furthermore, the 
adopted higher Ztmn are roughly consistent with the latest 
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Figure 1. Time evolution of spatial distributions of "FS" (field 
star) particles (cyan) and "GC" (globular cluster) ones (magenta) 
projected onto the x-z plane in the standard LMC model. For 
clarity, GC particles are plotted by bigger dots. The redshift (2) 
is given in the upper left corner for each panel. 



observation by WMAP suggesting that the epoch of reion- 
ization is 2: = IO.9I2 3 (Ps-gs 6t Bui. 2006). The models with 
Ztrun < 10 are found to yield too many GCs (> 1000) in the 
present models so that they can not be consistent with the 
observed number of old GCs (13) in the LMC. Therefore, 
we mainly discuss the results of the models with Ztrun = 15 
in the present study. 

It should be stressed here that without truncations of 
GC formation after reionization, the present models pro- 
duce too many GCs owing to a large number of low-mass 
halos virialized after reionization so that they can not be 
consistent with observations. This however does not neces- 
sarily mean that the truncation of GC formation by reion- 
ization is a crucial in better understanding the origin of the 
LMC's GCS: later, selective dynamical destruction of low- 
density GCs originating from low-density halos virialized af- 
ter reionization might well dramatically reduce the overpro- 
duced GCs in a model with no truncation of GC formation 
by reionization so that the simulated GC number can be 
consistent with the observed one even in such a model. Such 
later destruction of GCs is not modeled at all in the present 
study. 



2.4 Main points of analysis 

We investigate final structural and kinematical properties 
of stellar halos and GCs at 2; = 1 in models with different 
model parameters. We note that later accretion of satellites 
at 2 < 1.5 is minor in the final structures of the simulated 
halos at 2 = 0, so the calculation is ended at 2end = 1 
to obtain the dynamically relaxed halo structures. We first 
show the "standard" LMC model with Mt = 6.0 x lO^^M©, 
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5i = 0.39, A = 0.08, and Ztmn ~ 15 then discuss parameter 
dependences of the results. 

We adopt the initial LMC mass significantly larger than 
that of 2 X 10^" M0 used in the latest simulations for the 
orbital evolution of the present LMC interacting with the 
SMC (e.g., Yoshizawa & Noguchi 2003). The reason for this 
is that numerical simulations demonstrated that the mass 
of the LMC can become significantly smaller than its initial 
mass owing to stripping of dark matter halo and stars by 
the Galactic tidal field (Bckki & Chiba 2005): We need to 
adopt the LMC's mass significantly larger than the present 
LMC mass. The standard model shows that (1) the final 
mass within 7.5 kpc is about 2 x 10^" Mq and (2) the central 
velocity dispersion of the stellar halo is 50 — 60 km s^^ de- 
pending on projections. Models with smaller initial masses 
(~ 10^" Mq) shows that final velocity dispersions of stellar 
halos are significantly smaller than the observed one. 

Although the derived kinematics in the standard model 
is roughly consistent with observations (e.g., Minniti et al. 
2003), we investigate the following three models for com- 
parison: The "low density" model with Si = 0.19, the "high 
threshold" one with A'min = 1000, and the "low-z trunca- 
tion" one with ztrun = 10. The velocity dispersion of the 
stellar halo in the low density model (30 — 40 km s^^) is sig- 
nificantly smaller than the observed one, because the final 
single halo has a lower mean density and thus have a smaller 
meiss at z = 1 within 20 kpc that is more reasonable for the 
tidal radius in the more massive LMC at z = 1 than the 
present tidal radius of 15 kpc (van der Marel et al. 2002). 
Only subhalos with masses larger than lO^M© at the epochs 
of virialization can have FSs and GCs in the high threshold 
model. 

By assuming that the projected radial density profiles 
of GCSs (Sgc) are described as the power-law form; 

Egg oc R", (1) 

we derive the power-law slope (a) for each GCS. The derived 
power-law slopes can be compared with those of stellar ha- 
los with Sps oc R" so that we can discuss the differences 
in structural properties of these two old stellar components. 
The maximum velocities {Vm) and the central velocity dis- 
persion (do) in the radial dependences of rotational veloci- 
ties (Vrot) and velocity dispersions (cr) are derived to discuss 
kinematics of the two components. For convenience, Vm/crg 
is simply referred to as V/a throughout this paper. 

We adopt the slit size of 5 kpc (corresponding to the 
effective radius of the simulated GCS) so that we can esti- 
mate the radial profiles of Kot and cr with reasonably small 
error bars. Errors in Kot (o") are assumed to be equal to 
Kot/\/2(Af - 1) {cr/^2(N - 1)), where N is the total num- 
ber of particles for a given radial bin. Errors in V/a are esti- 
mated from the total number of particles at the radii where 
Kot becomes maximum. 

The total number of GCs are 70 - 450 (at z = 1) 
in the models with = 15 and thus much larger than 

the observed number (13) of old GCs in the present LMC. 
McLaughlin (1999) showed that total number of initial GCs 
in a galaxy can decrease by a factor of 25 within the Hubble 
time owing to GC destruction by the combination effect of 
galactic tidal fields and internal GC evolution (e.g., mass 
loss from massive and evolved stars). Therefore, it is reason- 
able to say that only several percent of the simulated GCs 



can survive to be observed as halo GCs in the LMC. We 
thus consider that the above range of GC number can be 
reasonable to be compared with observations. 

Table 1 summarizes the parameter values: Model name 

(column 1), Mt (2), A (4), Si (3), iV,nin (5), Ztrun (6), (^)pg 

(7), and (■^)qc (8). Here {^)^^ {{^)^^) is the mean V/a for 
the three projections {x-y, x-z, and y-z) for the stellar halo 
(GCS). Our previous simulations of disk galaxy formation 
shows that the final spin vectors of the simulated disks are 
similar to the initial spin vectors of dark matter halos (Bekki 
& Chiba 2001). We thus consider that the x-z plane of the 
simulated LMC corresponds to the disk plane of the LMC in 
order to discuss the observed kinematics of the stellar halo 
and the GCS. 

All the calculations have been carried out on the 
GRAPE board (Sugimoto et al. 1990). Total number of 
particles used in our simulations is 508686 and the gravi- 
tational softening length is 0.18 kpc for the LMC models. 
The adopted softening length is roughly similar to the ini- 
tial mean separation of the particles in a simulation. We 
used the COSMICS (Cosmological Initial Conditions and 
Microwave Anisotropy Codes) , which is a package of fortran 
programs for generating Gaussian random initial conditions 
for nonlinear structure formation simulations (Bertschinger 
1995). 

The present study does not intend to investigate de- 
struction of GCs and the resultant formation of field stars 
in the LMC's halo. The halo field stars originating from low- 
mass GCs and initially unbound or weakly bound star clus- 
ters (SCs) may well have structural and kinematical prop- 
erties similar to those of the GCS in the LMC, if the proba- 
bility of GCs (and SCs) being destroyed by the combination 
effects of the LMC's tidal field and internal GC evolution 
does not depend on orbits of GCs with respect to the LMC's 
center. As described later, the simulated halo and GCS in 
the LMC have similar dynamical properties. Thus, field halo 
formation via GC destruction would not change the main re- 
sults of the present models. 

As described later, the present models do not repro- 
duce self-consistently the observed kinematics of halo field 
stars and GCs in the LMC. This is in a striking contrast 
with previous simulations (e.g., Santos 2003; Bekki 2005) in 
which physical properties of the Galactic GCS can be well 
reproduced if truncation of GC formation by reionization is 
properly modeled. This is partly because the previous sim- 
ulations only tried to explain the Galactic GCS, where re- 
markable kinematical differences in halo field stars and GCs 
are not observed (i.e., they did not discuss stellar halos and 
GCSs in less luminous galaxies like the LMC). 



3 RESULTS 

3.1 The standard model 

Figure 1 shows the time evolution of spatial distributions 
of FSs and GCs from z = 30 to « = 1 that are formed 
within subhalos virialized before z = ztrun (=15). 434 small 
subhalos are virialized before z = ztrun and have masses less 
than 3.0 x lO^M© at the virialization and grow via hierarchi- 
cal merging with other subhalos with and without FSs and 
GCs {z = 8.6). These smaller subhalos with FSs and GCs 
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Figure 2. The final distributions of FS particles (upper) and 
GC ones (lower) projected onto the x-z plane in the standard 
model at 2 = 1. FS and GC particles form the stellar halo and 
the GCS, respectively, in the present study. The GCS appears to 
be slightly flattened in comparison with the stellar halo owing 
to the smaller number of GC particles with \z\ > 20 kpc. The 
overall distributions are however quite similar with each other in 
this model. 



>, 0.01 r 




R(kpc) 

Figure 3. The final projected radial profiles of FS (thin solid) 
and GC (thick solid) particles in the standard model. For com- 
parison, the profiles normalized to their central values arc shown. 
Both profiles can be well described as power-law ones with the 
slopes a ~ —2.0 for ii < 10 kpc. 



merge with one another to form bigger subhalos {z = 5.4), 
and finally these bigger halos also merge with one another 
{z = 3.2) to form a single halo till z — 1. FSs and GCs are 
tidally stripped from the subhalos during this hierarchical 
merging and consequently dispersed into the halo region to 
form a stellar halo and a GCS. 
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Figure 4. Radial dependences of rotational velocities Vrot (up- 
per) and velocity dispersions a (lower) for the stellar halo (thin 
solid) and the GCS (thick solid) projected onto the x-z plane in 
the standard model with the slit size of 5 kpc for the estimation 
of Viot a^iid ™ each bin. The projected distance here means the 
distance along the a;-axis in the simulation (See Figure 2). The 
results are shown for kpc ^ x ^ 25 kpc in the Vrot profile and 
for all selected particles with \x\ ^ 25 kpc in the a one. Viot = 
arc plotted with no error bars for bins with no GC particles (e.g., 
at a; ~ 25 kpc). Although the error bars are not small, it is clear 
that both the stellar halo and the GCS show little rotation. 
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Figure 5. The same as Figure 4 but for the halo and the GCS 
projected onto the y-z plane in the standard model. 
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Figure 2 shows that both the stellar halo (composed of 
FSs) and the GCS (GCs) have similar spherical distribu- 
tions in the x-z projection, though the GCSs has a smaller 
number of GCs with \z\ > 20 kpc. There are no significant 
differences in distributions projected onto the x-y, x-z, and 
y-z planes between the two components: both components 
show flattened distributions in the three projections with 
the major axes aligned with each other. These similarities 
means that the internal structures between the two com- 
ponents arc almost identical. The half-rmmber radius is 5.0 
kpc for the stellar halo and 5.3 kpc for the GCS, which sug- 
gests that there are no significant differences in dynamical 
properties between those two components. 

Figure 3 shows that both the stellar halo and the GCS 
have radial density profiles that can be approximated by the 
power- law ones with the slopes of a ~ —2 at least for i? < 20 
kpc. The apparent lack of flattening in the proflle of the stel- 
lar halo in the inner part of the simulated LMC {R < 2 kpc) 
is inconsistent with the best-fit exponential profile by Alves 
(2004), which shows flattening (or "core") of the proflle for 
R < 1 degree from the LMC's center. Although this incon- 
sistency between the simulated and the observed halos of 
the LMC has some profound physical meanings about the 
formation processes of the stellar halo (Bekki & Chiba 2007, 
in preparation), we intend to discuss this not in this paper 
but in our future papers. 

Figures 4 and 5 show that there are no remarkable dif- 
ferences in the radial dependences of rotational velocities 
(Kot) and velocity dispersions (a) between the stellar halo 
and the GCS. Both components show overall small Kot 
(< 40 km s"'^), radially decreasing a profiles, and small 
V/a (< 0.6). The estimated V/a in the x-y, x-z, and y- 
z projections arc 0.25±0.18, 0.57±0.13, and 0.21±0.01, re- 
spectively, for the stellar halo, and 0.47±0.33, 0.60±0.10, 
and 0.31±0.13, respectively, for the GCS. The derived small 
V/a clearly indicates that both components arc dynamically 
supported by velocity dispersion rather than by global ro- 
tation. The results shown in Figures 3, 4, and 5 thus show 
that there are no significant differences in dynamical prop- 
erties between the two components, though they originate 
from different parts of subhalos virializcd before z — Ztmn- 

The GCS is formed from subhalos with different masses 
and epochs of virialization so that the spatial distributions 
of GCs originating from different subhalos are different with 
one another. For example, GCs from subhalos with masses 
more than lO'^M© at their virialization epochs have a half- 
number radius of 3.5 kpc and thus a more compact spatial 
distribution in comparison with the GCS composed of all 
GCs. Such a more compact distribution can be seen in the 
stellar halo composed only of FSs originating from more 
massive subhalos. Those results imply that stellar compo- 
nents formed in more massive subhalos at high z are more 
likely to be the inner parts of galaxies at z = 0. 

3.2 Parameter dependences 

The dependences of structural and kinematical properties of 
stellar halos and GCSs on model parameters are summarized 
as follows. 

(i) low-density model: Figure 6 shows that the final pro- 
jected distributions of the stellar halo and the GCS in the 
low-density LMC model with So = 0.19 appears to be flat- 



m 



I 

isi 



I 

ISJ 




-20 20 

X-axis (kpc) 

Figure 6. The same as Figure 2 but for tlie low-density LMC 
model. Note that both the stellar halo and the GCS appear to be 
flattened. 
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Figure 7. The same as Figure 4 but for the low-density LMC 
model. Vrot = (cr = 0) are plotted with no error bars for bins 

with no GC (or halo) particles (e.g., a; 20 kpc and x ^ 25 kpc 
in the Vrot profile of GCs). 
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Figure 8. The dependence of the total number of subhalos (Afjj) 
with GCs, 2vir > ztrun, and Mh > Mtij in the standard model 
with different M^^ for Ztrun = 15. Mj^ is the mass of a halo when 
the halo is identified as a "virialized system" in the simulation. 
Note that ATh is very small for Mth > 3 x 10^ M©. 



tened in comparison with the standard model. We confirm 
that this flattening can be seen in the three projections (i.e., 
x-y, x-z, and y-z) and thus suggest that the derived flat- 
tening is due to the lower initial density of the LMC in 
this model. The power-law slope a in the projected den- 
sity distribution is —2.5 both for the stellar halo and for 
the GCS, which means that the radial density proflles are 
slightly steeper in this model than in the standard one. 

(ii) low-density model: Figure 7 shows that the two com- 
ponents in the low-density model have a small amount of 
rotation (Kot < 20 km s^^), which suggests that the two 
components are dynamically supported by velocity disper- 
sion. The mean V/a is 0.60 for the stellar halo and 0.56 
for the GCS, which are slightly higher than those in the 
standard model (See table 1). Flattened shapes and small 
V/a in this model suggest that the two components have 
anisotropic velocity dispersions. 

(iii) high Nmin model: the high threshold model with 
Nmin = 1000 has more compact distributions of the stellar 
halo and the GCS with the half-number radii of the two 
equal to 4.5 kpc and 3.4 kpc, respectively. This is due to the 
fact that more massive subhalos virialized before 

^trun can 

finally settle in the inner region of the LMC in this model. 
The GCS appears to be more flattened than the stellar halo, 
though V/a is not different between the two components. 

(iv) low-z truncation model: both the stellar halo and 
the GCS have less compact spatial distributions and shal- 
lower radial density profiles with a ~ —1.5 for the central 10 
kpc in the low-z truncation model. Both the stellar halo and 
the GCS have small V/a (0.19 and 0.21, respectively), which 
means that the two components are supported by velocity 
dispersion more strongly. A large number of GCs (~ 5000) 
is totally inconsistent with observations, even if later GC 
destruction by galactic tidal fields are considered. 

(v) Thus the present models all show very small V/a 
(0.2 — 0.6) of GCSs, which is much smaller than the ob- 
served value of ~ 2. Small V/a can be seen in models with 
high spin parameters (e.g., A — 0.12) and low initial masses 
(Mt = 2.0 X lO^^M©) and thus suggests that the present 
dissipationless models can not reproduce well the observed 
rotational kinematics of the LMC's GCS. The inner flat- 



tened halo derived only in the low-density model can be 
consistent with the observed one by Subramaniam (2006) 
and thus implies that the LMC could be formed from a low- 
density galaxy-scale fluctuation. 



4 DISCUSSIONS 

Wo have shown that (1) there are no significant kinematical 
differences in the simulated stellar halos and GCSs in the 
LMC models and (2) V/a of the GCSs can not be as high 
as observed (^ 2) in the GCS of the LMC. These failures to 
reproduce the observe kinematical differences between the 
two old stellar components (i.e., stellar halo and GCS) in 
the LMC imply that the adopted models for the formation 
of the two components at high redshifts can lack some im- 
portant ingredients of GC formation. These failures can be 
due to the adopted assumptions that (i) all subhalos with 
different masses and different redshifts of virialization can 
have both field stars and GCs and (ii) the GCS was formed 
from dissipationless merging of subhalos that had been viri- 
alized before reionization (i.e., Zvir > Ztmn) and thus had 
GCs (before their merging leading to the formation of the 
LMC at later redshifts). 

Observational studies of GCs for galaxies in the Local 
Group of galaxies showed that dwarf galaxies fainter than 
Mv = — 13 mag appear to have no GCs (van den Bergh 
2000). This result means that (i) there could be a possible 
threshold galaxy mass (10* - lO^Mo) above which GC for- 
mation is possible and (ii) the present dissipationless mod- 
els discussed so far did not consider this possible threshold 
mass (Mth)- We accordingly investigate how the total num- 
ber of subhalos that can be virialized before ztmn and thus 
have GCs (A'h) depend on Mth for the standard model with 
^trun = 15. Figure 8 shows that (i) iVh is smaller for the 
models with larger Mth and (ii) if Mth > 3 x lO^^M©, al- 
most no subhalos virialized before ^trun = 15 can have GCs. 
The required Mth for no halo formation with GCs before 
reionization is higher for lower Ztrun- 

These results in Figure 8 therefore imply that if there 
is a threshold halo mass (Mth) for GC formation and if 
Mth > 3 x IO^Mq, the LMC's GCS can not be formed from 
hierarchical merging of subhalos with z^ir > «trun and with 
GCs. As shown in the present studies, dissipationless merg- 
ing of subhalos is responsible for the larger velocity disper- 
sion of old stellar components. These results in Figure 8 ac- 
cordingly imply that the observed small velocity dispersion 
of the LMC's GCS is due to the fact that the LMC's GCS 
was not formed from dissipationless merging of low-mass 
subhalos with Zvir > Ztrun- Figure 8 thus implies that old, 
metal-poor GCs in the LMC were not formed in low-mass 
subhalos that were the building blocks of the LMC, because 
the masses of the subhalos were systematically lower than 
Mth before reionization. 

If the LMC's old GCs do not originate from subhalos 
virialized at high redshifts, how were they formed ? Previ- 
ous numerical simulations showed that GCs can be formed 
during dissipative merging between the Galaxy and gas-rich 
dwarfs (Bekki & Chiba 2002) . We accordingly consider that 
dissipative merging of gas-rich subhalos can trigger the for- 
mation of GCs at the very early epoch of the disk forma- 
tion of the LMC. Formation of GCs via dissipative merg- 



Kinematics of the LMC 9 



ing of subgalactic clumps (e.g., gas-rich dwarfs) in the very 
early epoch of galactic disk formation results in disky spatial 
distributions and rotational kinematics of GCSs (Bekki & 
Chiba 2002). If old, metal-poor GCs in the LMC are formed 
by the above dissipativo processes, they should have slightly 
younger ages than the Galactic counterparts. 

We thus suggest the following possible scenario for the 
origin of the observed kinematicaJ difference in the stellar 
halo and the GCS of the LMC: (1) the stellar halo was 
formed from merging of low-mass subhalos with field stars 
and with no/few GCs (i.e., GC-less galaxy building blocks) 
and thus shows a large velocity dispersion and a low V/a 
and (2) the GCS was formed through dissipativo merging 
of gas-rich subhalos and gas infall at the very early epoch 
of the disk formation and thus shows rotational kinematics. 
In the first part (1) of this scenario, the masses of subhalos 
virialized before reionization were well below the threshold 
mass (Mth) for GC formation so that the halos could not 
form GCs. We thus suggest that Mth can cause differences 
in structural and kinematical properties between stellar ha- 
los and old GCs in less luminous galaxies like the LMC. 

Olsen et al. (2004) found that the kinematics of GCSs 
in lato-typc galaxies in the Sculptor group are consistent 
with rotational kinematics seen in HI components of these 
galaxies and suggested that the GCSs were formed in disks 
rather than in halos. Beasley et al. (2006) also found a large 
V/cT ~ 3 of the GCS in the low-mass dwarf galaxy (VCC 
1087) in the Virgo cluster of galaxies. These observations 
imply that the rotational kinematics seen in the LMC's GCS 
is not exceptional but can be found in GCSs of many less 
luminous galaxies and thus that the formation processes of 
GCSs in these galaxies can be discussed in terms of the pro- 
posed scenario above. We also suggest that there can be two 
different formation processes of old, metal-poor GCs before 
and after reionization: one is GC formation in high-density 
central regions of subhalos early virialized before reioniza- 
tion and the other is GC formation in the very early stage of 
disk formation. Future observational studies on the shapes 
of the stellar halos in these galaxies with GCSs having rota- 
tional kinematics will enable us to discuss the origin of the 
possible differences in structures and kinematics between the 
stellar halos and the GCSs for these galaxies in the context 
of the proposed scenario. 



the adopted assumption that both the field stars and the 
GCs in the LMC were formed within all low-mass subhalos 
virialized before reionization is highly likely to be wrong. 

We accordingly considered that there could be a thresh- 
old halo mass (Mth) above which GCs can be formed (i.e., 
below which only field stars can be formed) and investigated 
how the number of subhalos (A^h) that can be virialized be- 
fore ztTun and have GCs depends on Mth for ztiun = 15. 
We found that if 3 x 10* Mq ^ Mth, the present LMC can 
not contain GCs formed within subhalos with the rcdshifts 
of their virialization («vir) larger than that of reionization 
(^^trun). We also suggested that if «trun is lower, the required 
Mth (for no GC formation) needs to be higher. We there- 
fore concluded that if old, metal-poor GCs in the LMC have 
rotational kinematics, they were not formed in subhalos viri- 
alized before reionization. 

We suggested a possible scenario in which the stellar 
halo was formed from low-mass subhalos virialized before 
reionization and having no/few GCs whereas the GCS was 
formed at the very early epoch of the disk formation via 
dissipative merging of gas-rich subhalos and gas infall well 
after reionization. In this scenario, the origin of the observed 
possible rotational kinematics of the LMC's GCS is closely 
associated with dissipativo gas dynamics in the disk forma- 
tion of the LMC. The LMC's GCs are thus suggested to be 
slightly younger than the Galactic counterparts. It is how- 
ever unclear how old GCs can be formed during dissipative 
formation of the main body of the LMC. We thus plan to 
investigate whether the observed possible global rotation of 
the GCS in the LMC can be reproduced by our more so- 
phistical numerical models with gas dynamics and GC for- 
mation. 

We also suggested that the threshold halo mass (Afth) 
for GC formation can cause significant diS'erences in struc- 
tural, kinematical, and chemical properties between stellar 
halos and GCSs in less-luminous galaxies like the LMC. Fu- 
ture observations will extensively investigate structural and 
kinematical differences in stellar halos and GCs for galaxies 
beyond the Local Group and thus confirm the presence or 
the absence of the differences. We plan to investigate depen- 
dences of physical properties of stellar halos and GCSs on 
physical conditions of their host galaxies at their formation 
epochs (e.g., masses and spin parameters) based on fully 
self-consistent chemodynamical simulations. 



5 CONCLUSIONS 

We numerically investigated structural and kinematical 
properties of the stellar halo and the GCS in the LMC 
by assuming that the two old components can be formed 
from dissipationless merging of subhalos that were virial- 
ized before reionization and contained both field stars and 
GCs. We particularly discussed whether or not the observed 
GCS's rotational kinematics {V/a ~ 2) can be reproduced 
by the present models. Our simulations with difi'erent model 
parameters showed that the GCS composed of metal-poor 
GCs formed high redshifts (z > 6) before reionization has 
little rotation {V/cr ~ 0.4) and structures and kinematics 
similar to those of the stellar halo. This inconsistency there- 
fore implies that if old, metal-poor GCs in the LMC have 
rotational kinematics, they are highly unlikely to originate 
from the low-mass subhalos that formed the stellar halo: 
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